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ABSTRACT 

A new solution method is presented for steady-state, momentum-conserving, 
non-axisymmetric bow shocks and colliding winds in the thin-shell limit. This is a 
generalization of previous formulations to include a density gradient in the pre-shock 
ambient medium, as well as anisotropy in the pre-shock wind. For cases where the 
wind is unaccelerated, the formalism yields exact, analytic solutions. 

Solutions are presented for two bow shock cases: (1) that due to a star moving 
supersonically with respect to an ambient medium with a density gradient perpendicular 
to the stellar velocity, and (2) that due to a star with a misaligned, axisymmetric wind 
moving in a uniform medium. It is also shown under quite general circumstances that 
the total rate of energy thermalization in the bow shock is independent of the details 
of the wind asymmetry, including the orientation of the non-axisymmetric driving 
wind, provided the wind is non-accelerating and point-symmetric. A typical feature of 
the solutions is that the region near the standoff point is tilted, so that the star does 
not lie along the bisector of a parabolic fit to the standoff region. The principal use of 
this work is to infer the origin of bow shock asymmetries, whether due to the wind or 
ambient medium, or both. 



Subject headings: hydrodynamics — ISM: bubbles — ISM: HII regions — shock waves 
— stars: mass loss 



1. Introduction 

Supersonic stellar winds shock the surrounding gas and drive expanding bubbles into the 
interstellar medium. These shocks provide an opportunity to probe the properties of both the 
driving stellar wind and the ambient medium. If the star is moving with respect to the interstellar 
gas, the bubble will be distorted into a cometary shape. When the stellar motion is supersonic, 
we refer to these as stellar wind bow shocks (Baranov, Krasnobaev & Kulikovskii 1971; Dyson 
1975). Since the discovery of such bow shocks around young B stars (Van Buren & McCray 1988), 
bow shocks have been found associated with many classes of objects, such as pulsars (Kulkarni 
et al. 1992) and cataclysmic variables (e.g. Vela X-l: Kaper et al. 1997); examples include 
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well-known naked-eye stars (e.g. Betelgeuse: Noriega-Crespo et al. 1997). Bow shocks have been 
proposed as an explanation for cometary, ultracompact HII regions (Van Buren et al. 1990; Mac 
Low et al. 1991) and as a means of explaining the lifetimes of ultracompact HII regions. In a 
recent survey of the IRAS database using HiRes processing, Van Buren, Noriega-Crespo, & Dgani 
(1995) found 58 candidate bow shocks. 

Non-axisymmetric stellar wind bow shocks occur when a star with an anisotropic wind moves 
supersonically with respect to the local medium, or if the star has an isotropic wind but moves in 
an ambient medium containing a transverse density gradient. Models of non-axisymmetric bow 
shocks are relevant to cometary ultracompact HII regions due to wind-blowing O stars moving 
supersonically with respect to the surrounding molecular cloud, when the ambient material 
does not have a constant density. A non-axisymmetric bow shock has also been invoked to 
explain the morphology of Kepler's supernova remnant, where the supernova ejecta collide with 
a non-axisymmetric bow shock generated by the pre-supernova wind (Bandiera 1987; Borkowski, 
Blondin & Sarazin 1992). Another example is the bow shock due to the head of a jet propagating 
into a region with a density gradient. Non-axisymmetric, ram-pressure balance models of the 
collision between a stellar wind and the photoevaporating flow from an externally illuminated 
circumstellar disk have been given by Henney et al. (1996). A formulation for steady-state 
non-axisymmetric bow shocks and colliding winds was given by Bandiera (1993). However, 
Bandiera's numerical method is sufficiently complicated that a simpler, analytic method is 
desirable. 

In this contribution, I present a method for solving the problem of steady-state, momentum- 
conserving, non-axisymmetric, thin-shell bow shocks and colliding winds. This is an extension of 
the previous analytic solution method of Wilkin (1996, hereafter Paper I) and of Canto, Raga, & 
Wilkin (1996, hereafter Paper II) to non-axisymmetric problems (see also Wilkin 1997a). 

An outline of the paper is as follows. In § 2, we formulate the problem of the steady-state 
collision of two winds, and in § 3 we treat the problem of a bow shock resulting from an isotropic 
wind interacting with a plane-parallel flow containing a transverse density gradient. In § 4, we 
allow for non-isotropic winds, especially an axisymmetric wind with random orientation of the 
symmetry axis with respect to the direction of stellar motion. The rate at which kinetic energy is 
thermalized for the bow shock is treated in § 5. Results and future directions of this research are 
summarized in § 6. 

2. Mathematical Formulation and Solution Method 

2.1. Description of the Collision Surface 

The hypersonic collision of two winds will in general result in a system of two shocks. A 
specific example, that of a stellar wind bow shock, is shown schematically in Figure |]. Because this 
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paper is concerned with steady-state solutions, the colliding winds are assumed to be unchanging 
in time. The stellar wind bow shock arises when a wind-blowing star moves supersonically with 
respect to the intersellar gas. In this case, we formulate the problem in the reference frame of 
the star, so the ambient medium is described as a wind of parallel streamlines impinging on the 
bow shock. For the collision of two winds in a binary star, we will neglect orbital motion in order 
to consider a steady-state problem in an inertial reference frame. In steady-state, the amount of 
mass and momentum within a given volume does not increase, and a flow pattern exists between 
the two shocks that carries away the mass and momentum deposited by the colliding winds. If 
the shocks are radiative, post-shock cooling lowers the temperature of the gas and leads to a 
large compression. We make the strong thin-shell assumption that cooling is so efficient that 
the shocked shell collapses to an infinitesimally thin layer, with a finite surface density a of 
matter. For this thin-shell assumption to apply, it is also necessary that the post-shock gas not be 
supported by magnetic fields, which, if well-coupled to the gas, could maintain a finite thickness 
even in the presence of efficient cooling. We assume there are no other forces such as those due to 
radiation or magnetic fields. The two shocked winds may in principle be separated by a contact 
discontinuity. However, there would then be supersonic shear across the interface that is expected 
to be unstable, leading to a mixing of the two fluids. A detailed treatment of the mixing is beyond 
the scope of this work. Instead, the mixing is assumed to be instantaneous, so the shell will 
have a unique velocity Vt at any location within it. This velocity represents an average of the 
turbulent fluctuations that would be present in a more detailed treatment. In steady-state, the 
geometric shape of the shell is unchanging in time, so the velocity of matter within the infinitely 
thin layer must be purely tangent to the shell. There will, however, be acceleration normal to 
the shell, because the fluid typically does not follow a straight path. The assumption that the 
incident streams are hypersonic, combined with the perfect cooling assumption, means that the 
flow within the shell will also be hypersonic, and pressure forces may be neglected in describing 
the motions along the shell.Q Defining a spherical coordinate system (r, 6,<j)), and denoting the 
radius of the surface by r = R(6,(f)), a complete description of this idealized shell is given by 
specifying the quantities R, a, Vt as a function of position (0, <p) within the shell. There are in 
fact only four independent quantities, because the condition that the motion within the shell be 
tangential implies that one need only solve for two velocity components. 



Clearly this assumption must break down for real, finite temperature systems near the stagnation point, because 
the tangential flow velocity in the shell vanishes at that point. However, pressure forces depend upon the gradient 
of the pressure, which also vanishes at the stagnation point, so the pressure forces may still be small compared to 
momentum deposition in the stagnation region (Wilkin 1997b). In any case, beyond a small region near the stagnation 
point, the tangential flow in the shell will be supersonic, provided the incident flows are supersonic and the shocked 
fluid cools efficiently. 
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2.2. Previous Treatment of the Problem 

The conservation laws of mass and momentum may be used to derive the properties of the 
shell in terms of those of the two incident winds. Since the momentum conservation law has three 
components, there will be four equations in four unknowns, or including the condition of zero 
normal velocity, five equations in five unknowns. Because the assumption of vanishing thickness 
eliminates one spatial variable, these equations will be partial differential equations (PDEs) in 
two spatial coordinates. Suitable boundary conditions must also be specified. In practice, this 
means identifying the stagnation (or standoff) point, where the two winds collide head-on, and 
where for steady-state, the ram pressures of the two fluids balance. The equations may then be 
integrated away from the standoff point, following the motion of a fluid element. In order to begin 
the numerical integrations, one generally finds that an expansion about the conditions at the 
stagnation point is needed. 

This approach was taken by Bandiera (1993), who derived a set of PDEs in curvilinear 
coordinates matched to the shape of the shell. These equations were then solved under the 
assumption of radial, constant velocity, isotropic winds from two point sources. The moving stellar 
wind bow shock problem is formally obtained by taking the limit of one source placed infinitely 
far away, while allowing its mass loss rate to be infinite so as to produce a finite density near the 
other star (at the bow shock). 

Bandiera noted that for the specific problems cited, the motion within the shell would be 
along planes. One may readily see this for the bow shock shown in Figure 1. Consider a point 
on the bow shock surface, and the plane of constant azimuthal angle containing it and the stellar 
trajectory. The radial wind striking the shell at this point has momentum lying in this plane. 
Similarly, the ambient medium striking this point has momentum lying in the plane. If the two 
incident streams mix instantaneously on impact, the resultant momentum must also lie in the 
same plane. As the fluid flows along the shell, it continues to incorporate momentum contributions 
lying in the same plane. The collision of radial winds emanating from two point sources gives the 
same result. Thus, given the perfect mixing asumption, the flow within the shell will lie along 
planes, a situation we refer to as meridional flow. If the shocked fluid does not mix, the fluid 
trajectories in the shell are not confined to a plane, and further work is needed to solve this more 
complicated problem]^] 

The ensuing geometric simplification allowed Bandiera to construct a two-dimensional 
grid, consisting of neighboring trajectories within the shell at a set of azimuthal angles. His 
computational method was to integrate the PDEs for mass, normal momentum, and tangential 
momentum, by marching along the trajectories within the shell as if the equations were ODEs. 



2 In the case of axisymmetric problems with shear, one may show that for divergent flow fields such as those 
encountered in bow shocks and colliding winds, the flow geometry and the fluxes of mass, momentum and angular 
momentum are the same as for the mixed case (Wilkin, Canto and Raga 2000). 
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By differencing neighboring solutions, numerical values of the cross-streamline derivatives were 
obtained and supplied to the ODE integrator. Provided the spacing of the two-dimensional grid is 
sufficiently fine, his method will yield accurate solutions. Bandiera postulated that the dependence 
of the equations on derivatives of quantities across streamlines was "fictitious," although he did 
not succeed in eliminating it from the equations. 

In this contribution, the solution method will be simplified to an integral approach that avoids 
the need for PDEs, solving purely algebraic equations. Additionally, the solution may conveniently 
be obtained in terms of ordinary Cartesian, cylindrical polar, or spherical polar coordinates, and 
does not require a coordinate system matching the shape of the shell. 

2.3. The Solution Method 

The solution method is based upon the observation that thin shells driven by hypersonic 
winds are momentum-conserving in the vector sense (Paper I). In order to conserve momentum 
in steady-state, the momentum flux in the shell must be the vector sum of the two incident 
momentum fluxes, integrated from the standoff point to the point of interest. Such an integration 
is performed over the area of the shell between two planes of constant azimuthal angle <f>, in the 
limit that their separation A<p is infinitesimal (Fig. 1). For axisymmetric flow with no rotational 
motion about the symmetry axis, no mass or momentum crosses these bounding planes, so the 
mass and momentum flowing within such a wedge must exactly equal the sum of the mass or 
momentum fluxes onto the external faces of the bounding surfaces (shocks). Given the known 
vector momentum flux, one may determine the shape of the shell, since the fluid must move in the 
direction of its momentum. While this description is complete, the mathematics was subsequently 
simplified with the inclusion of an additional, angular momentum flux integral (Paper II). In 
Papers I and II, only constant velocity, isotropic winds were considered. For such winds, the vector 
momentum flux incident onto the shocked shell is independent of the detailed shape of the shell, 
so it is possible to specify the flux of momentum onto the shell from each side analytically. The 
methods of Papers I and II will now be extended to non-axisymmetric bow shocks and colliding 
winds, provided the flow is meridional. For meridional flow, the bounding surfaces are again planes 
of constant <f>, allowing us to integrate the external fluxes and determine the internal fluxes within 
the shell. For anisotropic, radial, constant velocity winds, such as those considered by Bandiera, 
the momentum flux will depend only upon the number of streamlines intersected, and we will be 
able to obtain solutions analytically. 

Proceeding with the solution method, we note the location of the stagnation point, where 
the two streams collide head-on. For the two-wind collision, this point lies between the two stars, 
along the line connecting their centers. For the bow shock due to a wind-blowing star moving 
with respect to the ambient medium, it will lie on the stellar trajectory, in the direction of stellar 
motion. For either problem, the z-coordinate axis is chosen to contain the stagnation point at 
6 = 0, with the coordinate origin at the wind source (Fig. 1). In what follows, the wind from 
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the coordinate origin will be referred to as "the wind" , while the "second wind" may be either a 
radial wind from a second point source, or the ambient medium in the frame of the moving star, 
where the ambient velocity is V a = —V a z. However, we will show explicit formulas only for the 
bow shock, reserving detailed treatment of the binary star colliding winds for a future paper. One 
should bear in mind that the methods shown below will be applicable as well for that problem. 
Returning to the description of a thin slice of the shell bounded by planes of constant 4>, the mass, 
momentum, and angular momentum flux functions per unit azimuthal angle are defined by 

$ m = waV t seer], (1) 

* = $m V t t, (2) 

$j = Rx*. (3) 

Here w = Rsm8 is the cylindrical radius, and the unit vector at constant 4> tangent to the shell is 
t = 4> x h/\(p x n|, where n is the unit outward normal to the shell. Noting that the arc length 
element ds traced along the shell is given by (ds) 2 = (dR) 2 + R 2 (d9) 2 + R 2 sin 2 9 (d(j)) 2 , the arc 
length element traced along the shell at constant 6 is w sec rj def), where the angle rj is given by 

dR 

tan?? = — /Rsm9. (4) 

ocp 

For example, the flux of mass crossing a surface of constant polar angle 6, between the azimuthal 
angles (j) and (j) + d(f), is $> m d(fr. If Vq is the ^-component of velocity within the shell, and we write 
the cylindrical components of the momentum flux function as Q m and & z , then these components 
are related to the angular momentum flux by 

*j = * m RV e i, (5) 

$j = R(& w cobO- ^ z sm6)4>. (6) 

For meridional flow, the angular momentum within the shell will be in the (^-direction, so for 
the remainder of the discussion we will only need the magnitude of the angular momentum flux 
$j = $j • cf>. 

The mass, momentum, and angular momentum flux functions due to the incident winds are 
similarly defined. For the first wind, located at the coordinate origin, these are given in terms of 
the wind density p w and velocity V„ as 

$ mtW A0 = j j Pw (V w - n) dA, (7) 

A0 = J J Pw (V w -n)V w dA, (8) 
j, w A<p = J J p w (V w ■ n)(R x VJ dA. (9) 
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The area of integration is that between two planes of constant <fi, from the standoff point to polar 
angle 9, following the shape of the collision surface, which is determined along the course of the 
integrations. The sign of the normal direction is chosen so as to point away from the origin. The 
flux functions for the ambient medium are strictly analogous, except the unit vector normal to the 
collision surface would be in the reverse direction: 

$ m , a A</> = ~ J J Pa(Va ■ n) dA, (10) 

* a A</> = ~ J J Pa(V a • n)V a dA, (11) 

*j, a A0 = -/|p a (V o -n)(Rx V a )dA. (12) 

If R(9, 4>) is the radius of the shell, and its partial derivatives are Rq and R^, the surface area 
element is 

dA = ^R 2 + R 2 e + Rl esc 2 9 R sin 9 d6 d(f>. (13) 
The unit outward normal n to the surface is given by 

{R r R$ 9 R(j> esc 9 0) 
n = = L . (14) 

^R? + R 2 + R 2 esc 2 9 

Combining these, 

hdA = (Rr -R e 9- R<p esc 6 $) R sin 9 d9 d<f>. (15) 

For a radial wind, the partial derivatives with respect to 9 and 4> do not enter the expression for 
V m ■ rid A, which is simply 

V w -hdA = V w R 2 sin 9 d9 d<j). (16) 

Similarly, for the ambient medium V a = — V a z = — y a (rcos6* — 9sin9), so noting that 
dw/d9 = R cos 9 + Rg sin 9, we have 

-V a -hdA = V a wdwd4>. (17) 

This eliminates some of the partial derivatives that complicated Bandiera's very general treatment 
of the geometry. The rate at which conserved quantities are advected onto the shell does not 
depend upon the detailed description of the shell, provided we ensure the correct number of 
streamlines is counted. Using eq.(16), the resulting forms of eqs.(7-9) for the radial stellar wind 
are 

$m, w = / R 2 p w V w sm9'd9', (18) 
Jo 



/ R 2 Pw V 2 [w sin 9' + z cos 9'\ sin 9' d9' , (19) 
Jo 

®j, w = 0. (20) 
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Here the radial unit vector has been written in terms of its cylindrical polar components, and the 
integrations are performed at constant cf>. Because the stellar wind is radial, it imparts no angular 
momentum to the shell, and the wind angular momentum flux vanishes. Using eq.(17) to simplify 
eqs. (10-12) for the ambient medium, 

$m,a = V a p a zu'dzu', (21) 

Jo 

* a = -V*z / p a w'dw', (22) 
Jo 

$>j, a = V^ r Pa w l2 dw'. (23) 
Jo 

The mass conservation law for the pre-shock radial wind in steady-state is 

^(r 2 Pw V w ^j =0. (24) 

Consequently, if the wind is of mass loss rate M w and streamline-average speed V w , its properties 
are given by 

M u 
47r r' 2 



PwV w = 7Z ^f w (0A), (25) 



Pw V 2 = ^^g w (0^). (26) 
47T r z 

Here the dimensionless functions f w and g w are normalized to have unit streamline-average over 
47T steradians. Note that eq.(26) assumes the stellar wind to be coasting. Given this assumption, 
g w is independent of the radius. The wind mass and momentum fluxes onto the shell may now be 
written 

M r 9 

$m,w = -r 21 / f w (9',(f>) sm 9' d0', (27) 

47T JO 

/ g w (9', <j)) [w sin 9' + z cos 9'} sin 9' d9' . (28) 
Jo 



_ M w V, 

^11) — 



Air 

In order to obtain steady-state solutions, the ambient density must be independent of the 
z-coordinate, although its dependence on the remaining zu, <p coordinates may be arbitrary. Thus, 
we assume an ambient density of the form 

Pa = P0 fa{™,<t>), (29) 

where po is the value of the ambient density along the stellar trajectory and /(0, 4>) = 1. The final 
forms of the ambient flux functions are thus 



$m,a = P0 V a / fa dm'. (30) 

Jo 

*a = -K,zV (31) 



-9- 



$J,a=P0K 2 / fa^dw'. (32) 





The mass, momentum, and angular momentum flux functions for both the wind and ambient 
medium are clearly streamline integrals that do not depend upon the detailed shape of the shell - 
the integration is essentially one over the solid angle of the A0 wedge, as seen by the origin. This 
is a consequence of the fact that the pre-shock media conserve these quantities in steady-state. 

Suppose we have performed the integrations for the mass, momentum, and angular momentum 
onto the narrow slice of the shell, for both incident winds, according to eqs. (20,27,28,30-32). In 
steady-state, these quantities do not accumulate at any location within the shell, but are carried 
away by the flow within the layer. The conservation laws of mass, momentum, and angular 
momentum take the form 

$m = $m,w + $m,a, (33) 

* = * w + * a , (34) 

$J = <S>J,w + $J,a- (35) 

Now we may describe the properties of the shell explicitly in terms of those of the external winds. 
Equation (6) yields the shell radius 

<& j 

R = — . (36) 

$ ro cos 8 - & z sin 8 

This equation combines the formulation of Paper II, which was in terms of azimuthally integrated 
fluxes, with the original treatment of Paper I, using fluxes in an infinitesimally thin wedge. In 
eq.(36), each momentum flux function within the shell is to be evaluated as the sum of appropriate 
source terms, using the conservation laws (eqs. [34,35]). It is to be stressed that if the flow is 
not meridional, eq.(5) does not hold, and the solution method is considerably more complicated. 
However, eqs. (6,36) are still valid provided the angular momentum flux is replaced by its 
component. From eqs. (1,2), the tangential velocity of matter in the shell is 

V t i = */* m , (37) 

while the mass surface density in the shell is given by 

a = (^/n7|*|)cos7/. (38) 



2.4. The Applied Torque Method and an Example 

To provide a specific example of the solution method, as well as a reference solution to be 
compared to in the following, we consider the simple problem of a bow shock from an isotropic 
wind in a uniform ambient medium of density p a . In this case, f w = g w = f a = 1, and the flux 
functions are 

$ m ,™ = ^(l-cos0), (39) 
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& w = M^Ye. (# _ s i n Q cos 9) + i sin 2 6] , (40) 
®j, w = 0, (20) 

<S>m,a = \fJ 2 PaV a , (41) 

*« = -^V/^, (42), 

$ J,a = \^PaVl (43) 

Because the stellar wind has uniform speed, we have replaced V w with V w . One could immediately 
substitute these results into eq.(36) to determine the shell's shape. However, before proceeding 
to apply the formalism, note that the equation for the shape of the shell depends upon a specific 
combination of three flux functions in eq.(36). One may divide the equation into two parts, one 
part for each wind, according to 

T fc = ro[$ ro , fc cot 6 - $ 2 , fc ] - $ J)fc , (44) 

and k = w or a for the wind and ambient (or second wind) sources. Physically, T represents the 
applied torque necessary to compress the fan of streamlines between the stagnation point and 
position R(9,<fi) into a unidirectional beam possessing the same momentum flux (we shall call 
T the required torque). This is clearly seen by considering the radial wind from the coordinate 
origin. It has no angular momentum, but an equivalent beam containing the same momentum 
flux, located at a specific point R(6,<p) does indeed have angular momentum about the origin, 
because the beam will not be radial. The required torque for the isotropic wind follows from 
eqs. (20,40,44): 

T w = - Mw ^ W (l-6cot9). (45) 

Using eqs. (41-44), the torque necessary to compress the ambient streamlines to a unidirectional 
beam is 

T a = -m*p a Vl (46) 

For the ambient medium, the required torque expression simplifies because the ambient medium 
has no cylindrical radial momentum, and its z-momentum flux is related to its mass flux by 
eq.(31), so in general, 

T a = wV a $ m , a -<5>. la . (47) 

Similarly, the radial wind's required torque is simplified due to its lack of angular momentum 
about the origin. 

The solution surface is now defined by 



T w + T a = 0. 



(48) 
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The interpretation of eq.(48) is that each wind supplies the torque necessary to compress the other 
wind. These torques are equal and opposite only if the correct value of the radius (or cylindrical 
radius w) is used, conserving angular momentum as well as linear momentum. The utility of this 
approach is that if we vary the properties of one of the winds, holding the other constant, the 
shape of the shell is affected only by the change in the applied torque from the second wind. 

Now the applied torque expressions for the wind (eq.[45]) and the ambient medium (eq.[46]), 
when substituted into eq.(48), quickly recover the solution of Paper I, 



R(9) = R esc 0^3(1 -0 cot 0), (49) 
where Rq is the standoff radius: 



M W V W 



(50> 



3. Solutions for Bow Shock with Ambient Density Gradient 

In this section we assume the asymmetry of the bow shock to be due to a density gradient in 
the ambient medium, while the the stellar wind is assumed to be isotropic. The next section will 
treat bow shocks from anisotropic winds. 

The mass and momentum contributions of the stellar wind to the shell are given by 
eqs. (39,40,20), and its required torque is given by eq.(45). The ambient medium is assumed to 
have a plane-parallel density stratification, so in eq.(29), we replace f a (w, (ft) simply by f a {x), with 
f a is normalized so that / a (0) = 1, and where 

x = zucos(j). (51) 

The integrated mass and angular momentum flux onto the wedge from the ambient medium are 

<f>m,a = PoV a / f a (w' COS 4>) w' dm'. (52) 
J 

rzu 

$J,a = PoV* / f a (m' cos <f>) m' 2 dm'. (53) 
Jo 

A complete description of the bow shock is now obtained by adding the contributions of both the 
wind and ambient medium using the conservation laws (eqs. [33-35]). 

3.1. Solution for Exponential Density Law 

We first obtain the solution for an exponential distribution of ambient density. Let the density 
scale height be H. Define for brevity the coordinate y = x/H, so the density law is 

f a (x) = exp(-y), (54) 
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and the mass and angular momentum flux integrals are 

$ m , a = H 2 p V a sec 2 <j> [1 - exp(-y)(l + y)}, (55) 

$ J,a = H 3 P0 V* sec 3 [2 - exp(-y) (y 2 + 2y + 2)] . (56) 

The ambient momentum flux follows from eq.(31), while the required torque follows from eq.(47), 
which yields 



T a = H 3 p V* sec 3 



y - 2 + exp(-y)(y + 2) 



(57) 



Using this torque formula and that for the stellar wind given by eq.(45), the shell's shape is given 
by eq.(48), which yields upon simplification 

/ 2 sec 2 0y- 1 [y-2 + exp(-y)(y + 2)] = i(l-0cot0). (58) 

Here I = H/Rq is the density scale height in units of the standoff radius. Letting y ax be the value 
of y = x/H appropriate for the axisymmetric solution of eq.(49), 

y ax = cos (f) ^3(1-9 cot 9), (59) 



eq.(58) takes the simple form 

1 

y 



y-2 + exp(-y)(y + 2) 



(60) 



This formula may be solved numerically for y(9,<p), so letting y s denote the solution for y to 
eq.(60), the shell's shape is given by 

R{9,4>) = H sec csc y a (6,<f>). (61) 

The result is a family of bow shock solutions distinguished by the value of the nondimensional 
parameter I. By examining eq.(60) we may deduce several properties of the solution. The left-hand 
side of the equation is monotonically increasing as a function of y. The signs of y and y ax must be 
the same, since it is due to the cos 4> factor, so it follows that y increases monotonically, although 
in a very non-linear fashion, as a function of y ax . Because the parameter I does not enter eq.(60), 
we see that there is one universal solution for the problem, although this solution is scaled by the 
value of Rq and stretched (distorted) depending upon the value of H/Rq. The solution for the 
radius at all angles 9, 4> and for all values of I follows from Fig. ^. In particular, as y — > oo, the left 
hand side of eq.(60) approaches unity, implying a breakout angle for the shell when y ax = v6, 
which yields 

1-6*00 cot 6>oo = 2 1 2 sec 2 <p. (62) 

The opening angle 9^ depends upon cos</>, and applies only to cos<^ > 0. In the tail of the bow 
shock, for cos0 < 0, we have y — * — oo. In this case, y depends logarithmically upon y ax . The 
most interesting part of the solution, for \y\ < 1 may be described by an expansion in terms of 
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small y ax ■ Noting that for cos = the solution is identical to that of eq. (49) , we anticipate that 
an expansion for small cos (f> corresponds to an expansion in terms of small y ax about the standard 
bow shock solution. Letting R ax be the value of R(9) for that solution (given by eq.[49]), we obtain 



R ~ Rn 



1 _13_ 2 J_ 3 11843 4 

+ 4 Vax + 160 Vax + 24Q Vax+ 1075000 Vax + "". 



The behavior of the solution near the stagnation point is given by 



R — -Ri 



1 | COS( t> ft | ( l | 13cos 2 2 



(63) 



(64) 



Unlike the axisymmetric bow shock, there is a linear term in the behavior near the standoff point, 
so it is not describable as parabolic with the z coordinate axis as the axis of a parabola, except 
for the special angles <j) = ir/2 or 3tt/2. As a consistency check, note that as the scale height 
increases relative to the standoff radius (/ = H/Rq — ► oo), the solution reduces to the standard 
axisymmetric bow shock (Paper I). The lowest order effect near the standoff point is a tilt of 
the bow shock head, described by the linear term in 9. This means that although the wind and 
ambient streamlines meet head-on at the standoff point, they are not normal to the shell at this 
location, quite different from the axisymmetric case. This effect is due to the instantaneous mixing 
we have assumed. If shear is present in the shell, the two colliding flows could have separate 
stagnation points, which would be located where the incident stream is indeed normal to the shell. 
These bow shock solutions differ from the standard model (Baranov, Krasnobaev & Kulikovskii 
1971) in that the exponential ambient density distribution implies a finite total mass on one side 
of the bow shock. The solution then more closely resembles those of two-wind collisions (Paper II) 
where there is a finite opening angle for the bow shock tail. Similar problems include blast waves 
in finite mass media (Koo & McKee 1990) and wind breakout from a stratified medium (Canto 
1980; Borkowski, Blondin & Harrington 1997). For the high density region, cos0 < 0, the shell's 
asymptotic shape in the tail region is 

2/ 2 zu 
zr; msec 2 <p exp(— — cos</>). (65) 

7T H 

The shell continues to be more distorted as it expands due to the exp(— y) factor. Because the 
shell has a finite opening angle on one side, not all wind streamlines intersect the shell, but those 
with 9 > 9 00 {(f>) freely expand to infinity. Further consequences of the finite opening angle of the 
shell are discussed in 85. 



3.2. Solution for Linear or Polynomial Density Law 

Consider a stratification of the ambient medium according to 

f a (x) = 1 + aix + a 2 x 2 + .... 



(66) 
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Of course, one must ensure that this expression is non-negative over the domain of interest. The 
fluxes of mass and angular momentum onto the shell are 



± r 



J, a 



™ PoV a 



- + y x + T x + ... 

1 ai 02 2 

- + T x + y x +... 



(67) 



(68) 



The ambient momentum flux follows from eq.(31), while the required torque follows from eq.(47), 
giving the result 



1 01 3 2 6 n+l! 

1 + —x H a 2 x + ... + 

2 10 n + 3) 



a„x n + 



Substitution into eq.(48), we obtain upon simplification 
- 2r 3a 2 



l + yZZ7COS0 + 



2 2 
W COS 1 



+ ... 



3(1 -0cot0). 



(69) 



(70) 



Restricting the treatment to a linear density gradient, with only a\ non-zero, this equation is 
cubic in the variable w. One may solve analytically or numerically for the function w(6, (f>), so the 
solution surface is then given by R(9, 4>) = vj(9, (f>) esc 0. Inclusion of higher order terms in the 
density law simply increases the degree of the polynomial, but the solution is obtained in the same 
way. 



The behavior of the solution near the stagnation point is given by 



R_ 

Rq 



a\Ro cos < 



+ 



1 , 5 2 
5 + ( 32" a? 



20 



a 2 )i?o cos2 ' 



+ 



(71) 



As a consistency check of the solutions, for the case of a uniform ambient medium {a\, 02, ... = 0), 
eqs. (70,71) reduce to the standard solution of Paper I. As a further check, note that this 
result is consistent with the solution for an exponential mass distribution, if we choose 
the coefficients corresponding to the exponential distribution of the previous subsection, 
ai = —1/H,a,2 = 1/2H 2 , recovering eq.(64). 



4. Treatment of Anisotropic Winds 

4.1. The Axisymmetric Wind 

We now relax the assumption that the wind is isotropic to permit an axisymmetric wind, 
where the axis of symmetry is misaligned with the direction of stellar motion, or in the case of a 
two-wind collision in a binary star, where the axis does not point to the other star. Treatment of 
the accelerated wind will be deferred to a future contribution; in this paper the wind is assumed 
to be coasting. The coordinate axes are chosen so that the z-direction points in the direction of 
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stellar motion, or towards the second star. In terms of a spherical coordinate system, where the 
azimuthal angle (p is measured about the z-axis, motion of the shocked fluid in the shell is along 
planes of constant 0. We also define starred coordinates so that the z*-axis is the symmetry axis 
of the stellar wind. The two coordinate systems are related by 

sin cos 0* = sin 9 cos 4>, (72) 
sin 9* sin^* = sin 6 1 sin ^ cos A — cos 6* sin A, (73) 
cos#* = sin 9 sin (f> sin A + cos 9 cos A. (74) 

The wind mass and momentum flux densities depend on the polar angle 9*: 

PwV w = -—^ f w (9*), (75) 



Pw V* = -^g w (e*)- (76) 



M w V, 

The nondimensional functions f w and g w are normalized to have unit average value over Air 
steradians. 

The incident fluxes of mass, momentum, and angular momentum from the wind onto the shell 

are 

*m,«,(M) = ^-F w (e,<f>), (77) 

«MM) = ^^G w (0A), (78) 

*j l1u (M) = 0. (20) 
where the nondimensional functions F w and G w = G w ^ w + G w z z are given by 



7? 
1 w 



= / f w sm9'd9', (79) 
Jo 
r e 

= g w [w sin 9' + z cos 9'] sin 9'd9'. (80) 
Jo 



Here f w and g w have argument 0*(6', (p), so the integrations are performed using the transformation 
eqs. (72-74) to evaluate 0* (#,</>). We also define a nondimensional function T w associated with the 
required torque T w to compress the wind streamlines according to 

%j = — 7— -wT w . (81) 
Because the wind's angular momentum flux vanishes, eq.(44) implies that 

T w = G ww cot 9 — G WjZ . (82) 
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The functions f w (9*), g w {9*) may now be expanded in terms of powers of cos#*, 

oo 

fw(0*) =X> COS '0*> (83) 

i=0 
oo 

9w(0*) = E c * cos ^*- ( 84 ) 

i=0 

Defining for brevity 

p = sin sin A, (85) 

q = cos A, (86) 

the transformation given by eq.(74) is cos#* = psin.9 + qcosO. Denoting the trigonometric 
integrals by 

l jk = / sm J 9' cos k 9' d9', (87) 
Jo 

we evaluate F w and G w using the above series, and bringing p and q outside of the integrals, we 
finally have 

oo i /A 

F w = E **Ef H ? j L j ^-H-iJ' ( 88 ) 



=0 j=0 



OO J 



G w = °i Yl P l 3 9 3 \-)[^ h+i-jj + z h+i-j,i+j]. (89) 



=o j=0 



Here we have used the binomial coefficients to write the sums. The applied torque necessary to 
compress the wind streamlines to a thin shell, in nondimensional form, is now 

oo 

T w = 5>rW. (go) 

j=0 

where the responses due to the individual cos L 9 if terms are given by 

= E jTV ( l ) [cot 9 \ 2+i _ j:j - h +i -j, 1+j ] . (91) 

j=0 \ J J 



4.2. General Solution for Bow Shock Driven by an Axisymmetric, Misaligned Wind 

The wind description of the previous subsection may now be applied to the problem of a bow 
shock driven by an anisotropic wind. The wind is assumed to be driven by a star moving at speed 
V a in a medium of uniform density p a . We describe the bow shock's properties in the frame of the 
star. 
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The standoff radius, denned as the shell radius at 9 = 0, which corresponds to 6* = A, is given 

by 

= s (92) 

V ^p a v* 

We also define Rq to be the standoff radius for the equivalent isotropic wind 



*WS (93> 



It is important to recall for the results below that the standoff radius will be Ro^g w (X), where 

oo 

Su, (A) =J2°i COS ' A - ( 94 ) 
i=0 

By eqs. (46,48,82), the solution for the shell's radius is 



R(9, (f>) = R esc 0^-QT w , (95) 

where T w is given by eqs. (90,91). The total mass and momentum fluxes, including the contribution 
from the ambient medium, are 

<S> m = ^L [Fw + (96) 

47T LOL 

* = ™ + (G w , z - \w 2 ) z], (97) 

where w = zu/Rq, and a = V a /V w . 

To obtain a complete description of the bow shock's properties, we must specify the velocity 
of material in the shell and the mass surface density of matter. The velocity is given by the ratio 
of the mass and momentum fluxes: 



[2 G w ^ w + (2 G w>z - w 2 ) z] 
[2aF w + w 2 



v t i = v a l w > w z:^ i\~ — • (98) 



The surface density is given by eq.(38), which yields 



a = p a R = cos 77. (99) 

2wJ[AG 2 +(2G^ z -w 2 ) 2 } 



In order to evaluate a, we need to know the angle rj. Using eqs. (4,95), we have 

1 1 d T 

tan^-csc*--^, (100) 

where differentiation of eqs. (90,91) with respect to (f>, which enters only in the variable p, yields 

fit 00 /yr(') 

-5* = £>-5#- (101> 
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cot - j)p % J q J ( . ) [cot 9h+i-j,j - 



3=0 



\3 



(102) 



Note that the summation limits are restricted so that the terms with j = i vanish, including the 
Co term. 

In the region of the standoff point, an expansion to second order in 9 gives 

'I P 2 ^ 2 (A) 



4s«,(A) 



5 32 5 2 (A ) 



S(p 2 g^(X)-qg' w (X)) 



40 5 .(A) jj- W 

Primes on the function g w represent differentiation with respect to q = cos A. As was the case 
of asymmetry resulting from an ambient density gradient, there is a 9 1 term indicating a tilt of 
the standoff region relative to the direction of the stellar motion. The mass and momentum flux 
functions near the standoff point are 



m w v w n 

47T 



(104) 
(105) 



^(A)^ + |<i(A)z}# 3 . 
Equations (4,28,29) now give the tangential velocity in the shell and the mass per unit area as 

y,t-Hy. <> fe^*+£A^* , (106) 
3 af w (\)+g w (\) 

and 



3Rop a {a f w (\) + g w (X)Y 



(107) 



4 ^g w (\)(gl(\)+ P '(g> w (W 
One may readily see that these solutions reduce to the results of Paper I for an isotropic wind 
(fw = 9w = 1). 



4.3. Simple Solutions: Quadratic Dependence on cos6* 

The simplest non-trivial solution is for a linear dependence of mass and momentum fluxes 
upon cos#*. However, in astrophysical applications one is frequently concerned with stellar winds 
that have symmetry with respect to the equatorial plane, such as when the asymmetry arises 
due to rotation of the star. For winds that are symmetric with respect to 9 = tt/2, the functions 
fw{9*), g w (9*) will have expansions in terms of even powers of cos 9* only. Thus, we shall give the 
solution for a wind that depends upon cos 2 #*. This is sufficiently general to include as a subset 
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previous models of non-axisymmetric bow shocks (Bandiera 1993; Chen & Huang 1997). We 
assume the mass and momentum fluxes are described by 



fw = bo + h cos + b 2 cos 2 0*, 

9w = CO + Cl COS 0* + C2 cos 2 0*, 
where the normalization requires bo = (1 — 62/3) and Co = (1 — C2/3). 
The mass flux integral i 7 ^ is 



(108) 
(109) 



F„, = j& (1 - /x) + -± \p (9 - sin 9 cos 0) + g sin 2 9] 
+ y sin 2 0[(g 2 -p 2 )cos0 + 2pgsin0]J, 



(110) 



where b' = bo + 62 (2p 2 + q 2 )/3. The components of G TO are 



G w , ro = i{c o (20-sin20) + | 



p (8 - 9 cos + cos 30) + 4g sin 3 9 



+c 2 shr 



2pg sin 9 + (q 2 — p 2 ) cos 9 



(111) 



G 



(( , 2 ^Ijc o (l-cos20) + ! 



4p sin 3 6* + g (4 - 3 cos 9 - cos 39) 



+C2 



^ (49 - sin 40) + ^ - P ^ (2 + cos 29) sin 2 



(112) 



Here the coefficient c = Co + C2(3p 2 + (7 2 )/4 depends on and A. Using eqs.(82,lll,112) we obtain 

T™ = {-^(l-0cot0)-|(l-cos0)[<7 + ptan(^)] 



+ ^ (p 2 -g 2 )sin 2 0-;pg(20-sin20) 
8 [ 

Using eq.(95), with T w given by eq.(113), we obtain the shell's shape 

R = Ro esc 0^3(1 - 9 cot 0)(c o + j(3p 2 + g 2 )) 

Q 

+2ci(l - cos0)[g + ptan(-)] 

3 1 1/2 

+^[(g 2 -;p 2 ) sin 2 + ^(20 -sin 20)] I , 



(113) 



(114) 



where the explicit dependence upon stellar inclination and azimuthal angle is obtained using 
p = sin</>sinA and q = cos A. Examples are shown in Figures 4-6 for even parity winds with 
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C2 = 3, 1.5, and —1.5. Noting that Bandiera's notation for this problem was Ay = — 2c2/3, these 
correspond to his cases of Ay = —2, —1 and +1. 

In the neighborhood of the standoff point, the shell's shape is given by 



The previous formulation is not restricted to an axisymmetric wind, so one may obtain 
solutions for bow shocks driven by radial, non-axisymmetric winds, by defining the momentum 
flux from the wind in terms of spherical harmonics. Because the method is clear from the previous 
discussion of an axisymmetric wind, we will not give further details here. We note only that the 
response from each spherical harmonic contribution to the wind momentum flux may be obtained 
from the required torque formula. Then the shape of the shell will depend upon the square root of 
the sum of these responses. 



The energy thermalization rate of axisymmetric bow shocks and colliding winds has been 
considered by Wilkin, Canto & Raga (2000). The thermalization occurs in the shocks, post-shock 
(radiative) relaxation, and mixing. The maximum possible rate of thermalization is given by the 
total incident kinetic energy flux in the center of mass frame. Here we are concerned with the 
total rate of energy thermalized, over the entire shell, rather than the amount per unit area. By 
center of mass, we mean the center of mass of the matter deposited onto the shell per unit time. 
This maximum thermalization rate is for the case of complete mixing. For the bow shock in a 
uniform ambient medium, the center of mass frame corresponds to the ambient frame, because a 
formally infinite amount of ambient mass per unit time strikes the shell, while the stellar wind 
has finite mass loss rate. For the bow shock driven by a star with a radial wind, the total energy 
thermalization rate is given by 



R 2 /Rl « [(c + aq + c 2 q 2 ) + |(ci + 2c 2 q)Q 
+(8c + 5c iq + 2c 2 (3p 2 + q 2 ))9 2 /20]. 
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4.4. The Non- Axisymmetric Wind 



5. Energy Thermalization Rates 
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Here primes indicate velocities in the ambient frame, so = V w f + V a z, giving 

K, 2 = [v 2 + 2v w v aCOS e + v 2 }. 



(117) 



The velocity of the shell is its pattern speed and corresponds to the stellar velocity V* = V a z. 
Because we are assuming a non-accelerating wind, the integral depends only upon the streamlines 
and not the detailed shape of the shell, so we may perform the integration over a spherical surface. 



For the case of an isotropic wind from a star moving through a plane-parallel stratified ambient 
medium, the rate of energy thermalization is precisely the same as for a uniform medium (Wilkin, 
Canto, & Raga 1997), since it is the kinetic energy in the ambient frame that is thermalized: 



(118) 



This rate of energy thermalization is appropriate for an arbitrary ambient distribution of matter, 
provided only that it is truly infinite in mass, and that the bow shock is steady-state in the star's 
frame. The second requirement demands that the ambient mass distribution be independent of 
the z-coordinate, although its (m, <p) distribution is arbitrary. An exception, for example, is the 
case of an exponential stratification, for which case the total ambient mass flux may be finite (on 
one side of the bow shock). For this case, the center of mass frame is not equal to the ambient 
frame, and it is important to separately define the center of mass frame for each cf> slice, as it will 
depend upon azimuthal angle for the portion of the bow shock that contains finite opening angle. 
Because of the resulting opening angle, not all of the stellar wind kinetic energy (in the ambient 
frame) will be thermalized, for two reasons: (1) some streamlines miss the shell, and (2) there is 
a net flow of momentum along the shell even at infinite distance, whereas for the standard bow 
shock Vt — > in the ambient frame, far in the bow shock tail. 

We now consider the case of an anisotropic wind, to determine whether the energy 
thermalization rate depends on the orientation of the wind. However, the thermalization rate 
remains independent of the stratification of the ambient medium, subject to the caveat about 
finite-mass systems. For the remainder of this section, we allow the wind to be non-axisymmetric, 
so in equations (75,76) we use f w = f w (9*, </>*), g w = g w (0*, </>*)• When we refer to an axisymmetric 
wind, its axis of symmetry will be 6* = as before. 

The total kinetic energy loss rate of the wind is 



ott Jo Jo Jw 

The functions f w and g w are assumed to have unit average over 4n steradians. In terms of this 
energy loss rate, we define the mean square wind speed as 



We also define a total vector momentum flux for the wind, a quantity that is non-vanishing only if 
the wind is not point-symmetric with respect to the origin 
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<V%>= 2E W /M W . 



(120) 




(121) 



r2w 



JO 



/ & w (l7,(f))d(f). 



(122) 
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We now wish to calculate the kinetic energy deposition rate to the shell in the ambient frame. 
Letting a = V*/V w , the wind velocity (squared) in the ambient frame becomes 



V' 2 = V 2 



.2 

Hi 

f 2 
Jw 



9v 

fw 



+ 2a^co S e + a 2 



The wind density is given by 



4ttR 2 p v 



M w f, 



The total kinetic energy flux onto the shell in the ambient frame is 



E, 



therm 



M W V 2 
8vr 



/*27T /*7T 

JO Jo 



^ + 2ag w cos8 + a 2 f u 

Jw 



sin 8* d8* d<j)* 



(123) 
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(125) 



The normalization condition for / permits the last term to be easily integrated. The second term 
is V* times the z component of the wind momentum flux. In terms of the total wind kinetic energy 
flux and momentum flux, we have 



Etherm = \m w {< V 2 > +V 2 ) + P w ■ V, 



(126) 



We see that a sufficient condition for the total thermalization rate to be independent of 
the wind orientation is the vanishing of P z>w . Because we do not wish to invoke a special 
alignment of the wind orientation with the star's direction of motion, more typically this requires 
a point-symmetric wind such that ~P W = 0. As a special case, this can be more explicitly 
confirmed for the axisymmetric wind, assuming it is symmetric with respect to its midplane 
8* = it/2. Because we are considering the total thermalization rate, we may perform the solid 
angle integration in starred coordinates. In this case, we need only the transformation equation 
for cos# which is given by eq.(74). We assume the wind to be axisymmetric and symmetric with 
respect to 8* = ir/2, so it has an expansion in cos 6** with only even powers. This implies that the 
2a cos 8 term in the integral will give no contribution, by parity, because it yields only a sin 0* 
term, which vanishes in the azimuthal integration and terms with odd powers of cost?*, which 
vanish in the polar angle integration. The remaining contribution is then 



-^therm — 



M W V W 



flu 



+ a 2 f w (8*) 



sin 8* d8*, 
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= -M w (< v w > +V 2 ). 



(128) 



The total energy thermalization rate is independent of the orientation of the midplane-symmetric, 
axisymmetric wind. For the case of a non-point-symmetric wind, including an axisymmetric wind 
that doesn't possess mirror symmetry with respect to the equator, the total thermalization rate 
will depend upon the stellar orientation because of the contribution from the 2a cos 8 term. This 
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is the case of the example solution given in Section 4.3, for the wind with quadratic dependence 
on cos 6*. For that wind, the total vector momentum loss rate is 

PM W V W „ / inft \ 
w = — g — ciz*, (129) 

where z* points in the direction of the wind symmetry axis z*. The thermalization rate for the 
bow shock driven by this wind therefor has a fluctuation depending on orientation of the wind 
symmetry axis with respect to the stellar velocity of magnitude 

AE therm = MwV ™ V * Cl cos x ^30) 

If the quadratic term in the wind momentum vanishes (c2 = 0), then the absolute value of c\ may 
not exceed unity. The amount this changes the total thermalization rate then depends upon the 
parameter a = V*/V w . 

We must note that changing orientation of the point-symmetric wind did not change the total 
thermalization rate, but the spectrum of shock emission will clearly be different. For example the 
peak post-shock temperature will depend upon the incident normal component of velocity. The 
sum of thermalization by shock, relaxation, and mixing is independent of orientation, but the 
individual contributions will vary. This highlights the fact that the properties of bow shocks due 
to non-axisymmetric winds may be determined by measuring the bow shock shape, mass surface 
density (or column density), kinematics, and radiated energy. Although detailed shock calculations 
are beyond the scope of this work, a substantial amount of information should be obtainable from 
these global quantities and may be sufficient for sources where the data are sparse. 



6. Summary 

I have shown how to solve the problem of non-axisymmetric bow shocks and wind collisions 
with a simple formalism that requires only algebraic equations. Most often one considers constant 
wind speed for such problems, and in this case the method leads to exact, analytic (although 
possibly implicit) solutions. The availablility of simple analytic solutions makes it much easier 
to model observed sources and derive the properties of the driving winds. Among the principal 
applications of these solutions would be to determine the cause of the asymmetry in observed bow 
shocks, whether it be due to the ambient medium or an anisotropic wind, or both. 

Future improvements are needed to include non-radial and accelerating winds and shearing 
motions in the shell, in which case the fluid elements are not restricted to a plane. Also, 
non-axisymmetric bow shocks due to colliding winds in binary systems, including the orbital 
motion, require a formulation in a non-inertial frame. A forthcoming paper will describe how 
to solve the general problem of colliding winds from two stars, including both anisotropy and 
acceleration in the winds. 
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Fig. 1. — Thin Shell Model for a Stellar Wind Bow Shock. The top panel defines the spherical 
coordinate system and shows a thin wedge cut by two planes of constant (f>, while the bottom panel 
shows the wind, ambient, and tangential flows in the frame of the star. 

Fig. 2. — Solution of eq.(60): Scaled x-coordinate versus that of a bow shock in a uniform medium. 

Fig. 3. — Solution surfaces for I = 3 (top) and I = 1 (bottom). Stars mark the source of the wind. 
The ambient density increases towards the left of the page, and the direction of stellar motion is 
downwards. The bow shocks are viewed from an angle of 10° by rotating about the x-axis. Contours 
of constant 6 are at every 8° up to 120°, while contours of constant 4> are at every 10°. 

Fig. 4. — Bow Shock for C2 = 3, a polar wind. The wind is inclined with respect to the z-direction 
by angles (left column, from top to bottom) 0,10,30 degrees, and (right column, top to bottom) by 
50,70, and 90 degrees. 

Fig. 5. — Bow Shock for C2 = 1.5, a polar wind. The wind is inclined by angles 0,10,30,50,70, and 
90 degrees with respect to the z-direction as in Figure 4. 

Fig. 6. — Bow Shock for C2 = —1.5, an equatorial wind. The wind is inclined by angles 0,10,30,50,70, 
and 90 degrees with respect to the z-direction as in Figure 4. 
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